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Abstract 

Wc prove, under some generic assumptions, that the semiclassical 
spectrum modulo 0{h^) of a one dimensional pseudodifferential oper- 
ator completely determines the symplectic geometry of the underlying 
classical system. In particular, the spectrum determines the hamilto- 
nian dynamics of the principal symbol. 

1 Introduction 

In this article I would like to advocate an inverse spectral theory for pseu- 
dodifferential operators. What does this means ? One of the most famous 
inverse spectral problems, made fashionable by Kac's very entertaining arti- 
cle [11], with a mind-catching title "Can one hear the shape of a drum ?"(^), 
was about the Laplace operator on a bounded domain 0, C M". Frequencies 
u solutions to the eigenvalue problem 

—Au = v^u, n = on dVt 
2 

may be viewed as harmonics that can be heard when the interior of the 
"membrane" vibrates freely. The question was whether the knowledge 
of all frequencies completely determines (up to isometry, of course). As 
Kac mentioned, this question appears naturally in the context of quantum 
mechanics, for a particle trapped in a hard potential well. An important 
observation in this paper was the relevance of the Weyl law, which let us 
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find the volume (or area when n = 2) of J7 from the asymptotic behaviour 
of large eigenvalues. 

Counterexamples are now known : there are non-isometric shapes in 

that produce different frequencies [9]. Nonetheless, this fact should not 
let us think that the problem has become obsolete. As far as I know, the 
seemingly simple case of a convex, bounded domain G with analytic 
boundary is still open (although by the time this article is published, it 
might very well have been settled by Zelditch, see [22, 21]). 

Understanding this problem requires putting it in a wider perspective. A 
natural variant of Kac's problem is whether the spectrum of the Laplace op- 
erator Ag on a compact riemannian manifold (M, g) determines the metric g. 
Although, here again, counterexamples have been known for a long time [14], 
our understanding remains relatively poor. Recent works by Zelditch and 
Guillemin suggest that microlocal tools are quite relevant for all these ques- 
tions. This, in turn, is a hint that more general operators than the Laplacian 
could be dealt with similarly. 

From a quantum mechanical viewpoint, Kac's situation is quite extreme. 
A more natural setting would involve a particle 'trapped' by a smooth poten- 
tial well. No more boundary problems, but instead a Schrodinger operator 
on R" 

P = -^A + Vix). 

Of course now the potential function V should be recovered from the spec- 
trum of P. This inverse spectral problem has been studied a lot, but only 
very recently have microlocal tools similar to those used by Guillemin and 
Zelditch been applied to it [10, 5]. 

Here, I would like to shift again the initial problem one step further away. 
Instead of the Laplacian, or the Schrodinger operator, why not consider any 
(elliptic) differential operator, or even, since we're at it, any pseudodifferen- 
tial operator ? Of course, since there is no domain $7 anymore, no potential 
function V , the sensible question is what should we try to recover from the 
spectrum ? 

The inverse spectral problems I've mentioned here can all be understood 
as semiclassical limits. From a quantum object, the spectrum, one wants to 
recover classical observables such as the metric g, or the potential V . These 
quantities, in turn, fully determine the classical dynamics of the system. 
For general pseudodifferential operators, semiclassical analysis still shows 
the strong relationship between the classical dynamics and the quantum 
spectrum, so I believe that the most natural "object" that we should try 
and recover from the spectrum is precisely this classical dynamics. This, 
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precisely, amounts to determining the principal symbol of the operator. In 
fact, if we keep in mind Weyl's asymptotics, this sounds fairly natural, for it 
is well known that Weyl's asymptotics extend to arbitrary pseudo differential 
operators, provided that we compute phase space volumes defined by energy 
ranges given by the principal symbol [15]. 

As in the riemannian case, one should take into account a symmetry 
group acting on the classical data. For general pseudodifferential operators, 
there's only one available : the group of symplectomorphisms, acting on the 
phase space M. This is a much bigger group than the group of riemannian 
isometrics, in accordance with the fact that the space of principal symbols 
C°°(M) is much bigger than the space of riemannian metrics, or potential 
functions. 

2 The setting 

Since we aim at recovering the classical dynamics from the spectrum, we 
are going to work in the setting of semiclassical pseudodifferential opera- 
tors, which we recall here. Throughout this work, we only consider the 
one- dimensional theory. It would be very interesting to have higher dimen- 
sional results, but it is not expected that such precise results would persist. 
However, a reasonable challenge would be to undertake a similar study for 
the completely integrable case. 

The classes '^'^{m) of semiclassical pseudodifferential operators we use 
are standard. Let M = T*M = ^^(j,^)- Let d and m be real numbers. Let 
S'^{m) be the set of all families ^))te{o,i] of functions in C°°(M) such 
that 

^C7„;i'^(i + |x|2 + |e|2)f, (1) 

for some constant Cq, > 0, uniformly in h. Then ^'^(m) is the set of all 
(unbounded) linear operators P on L^(M) that are h-Weyl quantisations of 
symbols p G S'^{m) : 

{Pu){x) = {Opl{p)u){x) = ei^--y'^^p{^,^;h)u{y) |«| . 

27rn J^2 

The number d in (1) is called the ?i-order of the operator. Unless specified, 
it will always be zero here. In this work all symbols are assumed to admit 
a "classical" asymptotic expansions in integral powers of h (that is to say, 
in the ladder {S'^{m))dez,d^do for some do S Z). The leading term in this 
expansion is called the principal symbol of the operator. 
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Thus, the Schrodinger operator P = — ^A + F on M is a good candidate, 
of /i-order zero, whenever V has at most a polynomial growth. 

We use in this article the standard properties of such pseudodifferen- 
tial operators. In particular the composition sends ^''^(m) x "^'^ {m') to 
^1'^+'^' {m + m'). Moreover ah P € *°(0) are bounded: L2(]R) ^ L2(]R), 
uniformly for < ^ 1. 

An operator P G ^{m) is said to be elliptic at infinity if there exists a 
constant C > such that the principal symbol p satisfies 

for |x|V|e|^ > C. 

If P has a real-valued Weyl symbol, then it is a symmetric operator 
on L? with domain C^(M). If furthermore the principal symbol is elliptic 
at infinity, then P is essentially selfadjoint (see for instance [7, proposition 
8.5])._ 

Finally, when P G ^'^(m) is selfadjoint and elliptic at infinity, then for 
any / E Cq°(M), the operator f{P) defined by functional calculus satisfies 
f{P) e nfcgN(*°(-^"i)). See for instance [7] or [16] for details. 

The advantage of the semiclassical theory is that it allows us to use 
richer versions of Weyl's asymptotics. Instead of considering the limit of 
large eigenvalues, we fix a bounded spectral window I = [Eq,Ei] C M and 
study the asymptotics of all eigenvalues in I, as h ^ 0. 

Definition 2.1 We say that Assumption A{P,J,I) holds whenever 

1. P is a selfadjoint pseudodifferential operator in "^^(m) with principal 
symbol p, elliptic at infinity. 

2. J d [0, 1] is an infinite subset with zero as an accumulation point. 

3. There exists a neighbourhood J of I such that p~^{J) is compact in 
M. 

If Assumption A{P, J, I) holds, we denote by T,fi{P, I) the spectrum of 
P = P{h) in / (including multiplicities). We denote by T,(P,J',I) the 
family of all {T,fi{P, I); h ^ J}. It is well known that T,fi{P, I) is discrete 
for h small enough (see eg. [16, Theoreme 3.13], in a slightly different setup). 
Notice that when m > 0, the properness condition 3. is always satisfied. 

Proposition 2.2 Let P be a selfadjoint pseudodifferential operator in ^(m), 
with principal symbol p, elliptic at infinity. Let J C M be a closed interval 
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such that p ^{J) is compact. Then for any open interval I C J there exists 
ho > such that the spectrum of P in I is discrete for h ^ ho. 

Proof. The case m > is probably the most standard. We recall it quickly. 

Case m > 0. — Let x £ Cq°{J) be equal to 1 on I. Then by pseudodif- 
ferential functional calculus, f{P) is compact for h small enough. Therefore, 
denoting by 11/ the spectral projector on /, we have that 11/ = n//(P) is 
compact. This implies that 11/ has finite rank : the spectrum in / is discrete. 

Case m ^ 0. — First we show that one can replace J by an unbounded 
interval containing /. Thus assume J is compact. For notational conve- 
nience we let J = [—1,0]. For any a g] — 1,0[, Sard's theorem ensures the 
existence of a regular value A € ]a, 0[ for p. Then C := p~^{X) is a compact 
1-dimensional submanifold of : it is a finite union of circles. Let Q be 
the unbounded component of \ C. Suppose first that > A. Since the 
differential of p does not vanish on C, p < A in all the bounded components 
of \ C. Therefore oo. A]) is compact and one may replace J by 

]— oo,A]. Now if on the contrary p^^ < A on the unbounded components, 
we have to apply the same argument for A' €] — !,«[. Because A' < A, the 
new bounded components contain the old one, and therefore consist of the 
points where p > A'. Then p~^([A', +oo[) is compact, and one may replace 
J by [A',+cx)[. 

For the rest of the proof, we may suppose that J =]— 00,0]. Let — e G 
J \ /, close to the origin. Let x S C°°(M) be such that 



Let pfi be the Weyl symbol of P and define pn. = X°Ph- Then pfi is a symbol in 
^{m), with pfi ^ — e/2 on and pn = pn outside the set p^ {]—oo, — e/3]). 
Because m ^ 0, the set ]— oo, — e/3]) is included in p~^{J) for h small 
enough and hence must be compact as well. Denote by P the Weyl quanti- 
sation of p. Then for h small enough (P — A) is invertible for all As/. We 
can write 



Since P — P is the Weyl quantisation of a compactly supported symbol (with 
support in a compact independent of h), it is of the trace class. By analytic 



'x{x) = -e/2 
< xix) = X 
^Xix) ^ -e/2 



for X G ] — oo, — e] 
for X ^ —e/3 
everywhere. 
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Fredholm theory, we may take the determinant of (Id+(P— A)~^(P— P)) and 
conclude that the spectrum of P consists of the zeroes (with multiphcities) 
of a non-vanishing holomorphic function and hence is discrete. □ 

The goal of this article is to recover the dynamics of the hamiltonian p 
in the region p"^{I) for any operator P for which Assumption A{P,J,I) 
holds, for some subset ^7 C [0, 1]. Of course if we can do it for an arbitrary 
compact interval / C M, we recover the full dynamics of p. 

It turns out that, under some genericity conditions, these inverse spec- 
tral problems are fairly easy, compared to the general multi-dimensional 
problems alluded to in the introduction, in the sense that they only require 
a few terms in the asymptotics of the spectrum. Having this in mind, for 
a S M we denote by J^, I) + 0{h") the equivalence class of all T,fi{P, I) 
modulo ff. Our main result is Theorem 5.2, but we also state several inter- 
mediate results that require weaker hypothesis. An informal statement of 
Theorem 5.2 is as follows. 

Theorem 2.3 (Theorem 5.2) Let Assumption A{P, J, I) hold, and de- 
note M = p~^{I). Suppose that p^M is a Morse function. Assume that the 
graphs of the periods of all trajectories of the hamiltonian flow defined by 
P\M, as functions of the energy, intersect generically. 

Then the knowledge of'E{P,J',I) -\- 0{h'^) determines the dynamics of 
the hamiltonian system p^M ■ 

In fact, we determine completely the Hamiltonian p up to symplectic equiva- 
lence. Perhaps the most difficult step, for which Weyl's asymptotics are not 
enough, is the seemingly simple question to count the number of connected 
components of p~^{E), for a regular energy E G I (Theorem 4.2). 

Although we state everything for pseudodifferential operators defined on 
on R, it is most probable that all results extend to the case of pseudodifferen- 
tial operators defined on a one-dimensional compact manifold equipped with 
a smooth density, and to the case of Toeplitz operators on two-dimensional 
symplectic manifolds. 

The plan of the paper follows a fairly logical progression. Since we always 
work modulo symplectomorphisms, it is not reasonable to look for a formula 
that would give the principal symbol p. Instead we will try to recover as 
many symplectic invariants as possible from the spectrum so that, given 
two spectra, we should be able to tell whether they come from isomorphic 
systems. 
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Thus, the geometric object under study is a proper map p : M — > M, 
where M is a symplectic 2-manifold. The simplest symplectic invariants of 
this map are in fact topological invariants, and are dealt with in Sections 3 
and 4. Indeed, it follows from the action-angle theorem that as soon as 
G M is a regular value of p, then the fibres of p consist of a finite number 
of closed loops, each one diffeomorphic to a circle. Therefore, we need to be 
able to detect 

1. Whether an energy iiJ G R is a regular or critical value of p; this is 
done in Section 3 (Theorem 3.1). 

2. When is a regular value, the number of connected components of 
the fibre p~"'^(£'); Section 4.1 discusses this point (Theorem 4.2). 

Putting these results together we are able to recover the topological type 
of the singular fibration (Theorem 4.5). Then in Section 5, relying on the 
classification result of Dufour-Molino-Toulet [8, 17] (and some additional 
argument) we finally manage to recover the symplectic geometry of the 
system (Theorem 5.2). 

3 Singularities 

In order to detect whether a given energy £"0 € M is a critical value of p or 
not, it is enough to know the spectrum of -P in a small ball around £"0, at 
least under some nondegeneracy conditions. 

Recall that a function / : M ^ M is said to have a nondegenerate critical 
point m € M when df{m) = and the Hessian f"{m) is a nondegenerate 
quadratic form. Since M has dimension 2, there are only two cases : 

1. Elliptic case : there are local symplectic coordinates (a;,,^) in T^M 
such that /"(?Ti)(x,^) = C{x'^ +C^), for some constant C 7^ 0. 

2. Hyperbolic case : there are local symplectic coordinates (x,^) in 
TjnM such that f" {rn){x,S) = Cx^, for some constant C 7^ 0. 

We refer to each of these two cases as the type of the singularity m. 

Theorem 3.1 Let I he an interval containing Eq in its interior, and let 
Assumption A{P, J , I) hold. Assume also that p has only nondegenerate 
critical values in I, and that any two critical points with the same singu- 
larity type cannot have the same image by p. Then from the knowledge of 
S(P, J, I) + 0(/i2) one can infer 
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1. whether Eq is a critical value of p or not; 

2. in case Eq is a critical value, the type of the singularity. 

This theorem is a corohary of the following proposition, where we consider 
the density of states in small regions around Eq. We could, equivalently, 
invoke Weyl's asymptotics. 

Proposition 3.2 Let 7 G (0, 1) and, for E & I, 

pn{E) = h^-^#{T.n{P,B{E,h^))). 
Then for any E ^ I, the limit p{E) = lim;i_>o Ph{E) exists (in [0, +00] and 

1. if E is a regular value of p, then p is smooth at E; 

2. if E is an elliptic critical value of p, then p is discontinuous; 

3. if E is a hyperbolic critical value of p, then p{E) = +00. 

Proof. In case E is a regular value, the result follows directly from Weyl's 
asymptotics, which in turn can be derived from a semiclassical trace formula 
as in [6], or from the semiclassical Bohr-Sommerfeld rules as in [18]. Let us 
recall the Bohr-Sommerfeld approach. There exists an e > such that the 
eigenvalues of P inside [E — €,E + e] modulo 0{h°°) are the union (with 
multiplicities) of a finite number of spectra Uk, k = 1, . . . ,N , where is the 
number of connected components of p~^{E), and each ak is determined by 
quasimodes microlocalised on the corresponding component. Precisely, the 
elements of ak are given by the solutions A to the equation 

g'^''\X;h) elnhZ, (2) 

where the function g^'^^ admits an asymptotic expansion of the form 

{X;h)^ gi'^ (A) + hgP (A) + h^g^'^ (A) + • • • (3) 

with smooth coefficients gj. Moreover, if we denote by Cfc(A) the A:-ieth 
connected component of p~^(A), in such a way that the family (Cfc(A)) is 

(k) 

smooth in the variable A, then g^ is the action integral : 

g^^\\)=l idx. (4) 

JCk(X) 



From (2) is follows that, for h small enough 

# K n B{E, e)) = (27rn)-i g^''\E + e;h)- g^'XE - e; h) 



+ 6, 
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where the (5 G [— 1 , 1] is here to take care of the appropriate integer part of 
the right-hand-side. Hence 



#Kn5(^,e)) = i2TTh)- 



2e 



dE 



+ 6. 



With e = , this gives 



#{akr\B{E,fif)) 



TT 



dE 



Summing up all contributions for k - 
theorem, with 

p{E) = 1 
vr 



1, . . . , A^, we get the first claim of the 



dE 



The second claim can be proved in a similar way, using Bohr-Sommerfeld 
rules for elliptic singularities [20]. For our purposes, a Birkhoff normal form 
as in [2] would even be enough, since we deal with energy intervals of size 
0{fV). Here again there exists an e > such that the eigenvalues of P inside 
[E — E + e] modulo 0{h°°) are the union (with multiplicities) of a finite 
number of spectra cj^. corresponding to the various connected components 
of p~^{E). The difference is that not all components need have critical 
points. In fact by assumption only one component may have an elliptic 
critical point. Let us call ak the corresponding spectrum, and Cfc(A) the 
corresponding family of connected components. Since an elliptic critical 
point is a local extremum for p, the sets Cfc(A) are empty for all A in one of 
the halves of the interval [E — E + e]. Without loss of generality, one can 
assume that Cfc(A) = 0, VA € — £,£'[. Then Ck{E) is just a point, while 
Cfc(A) is a circle for all A G ]E^ E + e]. 

The Bohr-Sommerfeld rules for elliptic singularities say that the elements 
of (Tfc are the solutions A to an equation of the form 



(5) 



where the function e^^^ admits an asymptotic expansion exactly as g^^^ 
above (3). What's more, it is equally true that the principal term is an 
action integral : 



e(^)(S) = 0, 



(A) 



idx, \IX£[E,E + e]. 
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□ 

Calculating along the same lines as above, we find, for the quantity 

/,jf)(A) := ;ii-^#KnB(A,r))), 

the following limits : 

1. when \&[E- e,E[, limn^o Pn\^) = 0; 

2. when X€]E,E + e,E], lim^^o ( A) = ^ 

Finally, let £" be a hyperbolic critical value for p. Weyl asymptotics 
for such a situation have been worked out in [1], and the singular Bohr- 
Sommerfeld rules have been established in [3]. Using the latter result it can 
be proven as in [12] that the number of semiclassical eigenvalues generated 
by a hyperbolic fixed point, in a neighbourhood of size e = of the critical 
value, is of order e \\nh\ /h. Therefore, since there may be only one hyper- 
bolic point in p~^(0), it follows from this estimate and the results we just 
proved above for the regular and the elliptic case that 

pniE) ^C\lnh\, 

for some constant C > 0. This gives p{E) = +oo. 

Remark 3.3 It is probable that the nondegeneracy condition can be 
avoided. It is known quite generally that Weyl asymptotics hold for critical 
energies [23]. Thus, in all case, we recover the action integral as the inte- 
grated density of states. It would remain to show that the behaviour of the 
action integral completely determines the singularities of p. This is easy in 
the Schrodinger case p = C'^ + V{x). A 

4 Topology 

As we already mentioned above, once the singular fibres of p have been 
excluded, the topology is easy to understand. The map p become a locally 
trivial fibration whose fibres are disjoint unions of circles. 

Thus, if £^0 is a regular value of p, the semiglobal problem around Eq 
just amounts to counting the number of connected components of p~^{Eq). 



dX ' 



dE ' 
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The topology of singular fibres strongly depends on the type of singu- 
larity. Under the nondegeneracy assumption, the topology of the singular 
foliation in a neighbourhood of a singular fibre is essentially determined by 
the type of the singularity, and thus by Theorem 3.1. 

4.1 Connected components 

Let / be a compact interval of regular values of p. As above, we denote 
by Cfc(A), for k = 1,...,A^ and A G / the smooth families of connected 
components of p~^(A). Each Cfc(A) is globally invariant by the hamiltonian 
flow generated by p. Thus, this flow is periodic on Cfe(A). Let |Tfc(A)| 7^ 
be its primitive period (the sign is determined by the formula below). It 
follows from the action- angle theorem that is a smooth function of A. In 
fact it is well know that the period is the derivative of the action, and we 
have already met this quantity in the proof of Proposition 3.2. Using the 
action integral (4), we get 



Notice again that Tk never vanishes on /. 

Definition 4.1 We say that a point (A, t) ^ {I y. R*) is resonant whenever 
there exist {k,j) and {k',j') in {1, . . . ,N} x Z*, with k 7^ k' , such that 



Theorem 4.2 Let I be an interval of regular values of p, and let Assump- 
tion A{P, J, I) hold. Assume also that the set of resonant points in I xM 
is discrete. Then the number N of connected components ofp~^{\), As/, 
is determined by spectrum JT", /) + 0{h?). 

Before proving the theorem, let us just remark that the leading term of 
Weyl's asymptotics is not sharp enough for this. Indeed, it only gives the 
density p (Proposition 3.2) : 



From this one cannot distinguish, for example, one component with period 
T from two components with periods |ri| + \t2\ = \t\. 



rfc(A) 




dX 



jTkiX) =j'Tk'{X) = -t. 




(6) 



k=i 
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Remark also that the condition on resonant points is not adapted to 
systems with symmetries. For instance, a Schrodinger operator with a sym- 
metric double well has two components with equal periods. 

Proof of Theorem 4.2. We introduce the period lattice Ck{I) ■ 

Ck{I) ■ = {(A,i) G / X M; exp(i^p) is periodic on Cfc(A)} 
= {(A,jrfc(A)); A G 7,i G Z} 

and £(/) = U^i ^k{I)- The set C{I) is a union of smooth graphs that may 
intersect. The intersection points for t ^ Q are precisely the resonant points. 

In order to prove the theorem, we split the argument into two steps. The 
first one is to prove that J:{P,J,I) + 0(f?) determines £(/). The second 
step consists in showing why the knowledge of >C(/) — and the assumption on 
the set of resonant points — allows us to count the number of connected 
components. 

Step 1. Coming back to the Bohr-Sommerfeld rules discussed in the proof 
of Proposition 3.2, we recall that the spectrum of P modulo 0{h°°) is the 
superposition of the spectra generated by C^, for k = 2,...,N. For 
each k, Ck has a periodic structure that makes it close to an arithmetic 
progression. Thus, a simple and naive idea to distinguish between the dif- 
ferent periodic structures is to perform a frequency analysis, via a Fourier 
transform. Because we have at our disposal only a truncated sequence of 
eigenvalues (those that belong to /), we need to introduce a cut-off. Let 
I' (£ I and let x ^ C°°(M) have compact support in the interior of / and be 
equal to 1 on /'. We introduce the spectral measure 

Do{X;h)= Yl x{E)5e{X), 

where 5e is the Dirac distribution at E. The quantity we want to investi- 
gate is its Fourier transform. Since the mean spacing between consecutive 
eigenvalues is of order h, we use a corresponding scale for the time variable 
t, and thus introduce 

Z{t;h)= Yl 

The function Z is called the partition function. In fact, the idea we've just 
described is very well known in the semiclassical context, and is part of 
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the general formalism of trace formulce. We can consider the Schrodinger 
group U{t-h) = exp{-itP/h), and then Z{t;h) = Trace{x{P)U{t; h)). It is 
well known that xiP)U{t; h) is a Fomier Integral Operator, whose canonical 
transformation is the classical flow of p. Moreover, its trace is a lagrangian 
(or WKB) distribution associated with the lagrangian manifold of periods 



Such a result would almost finish the proof of Step 1. In fact, this state- 
ment exists in many versions, depending on various possible situations and 
hypothesis. For this reason we are not using it here as is, but instead resort 
once again to the Bohr-Sommerfeld rules, which is arguably the easiest way 
to go. 

We can split the partition function as 



Ap = {{E,t) G M^; 3z ep-\E),exp{TXp){z) = z} = C{I'). 



N 



k = l EGCTk 



Then from (2) one can introduce c i-^ f^^\c;h) as the inverse of A 
g^'^^X; h), which exists for h small enough, and write 




(7) 



(which, as before, is a finite sum) with 



N 




k=l 



Note that ipt{-;h) € C^(M). By the Poisson summation formula. 




(9) 



(which, contrary to (7), is a truly infinite sum) with 




(10) 



and 
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The integral is a compactly supported oscillatory integral, whose phase 
is stationary when j + 1 = or equivalently 

t = -m{\), x = f^''\c;h). (11) 



Moreover, the Hessian of the phase, Q := t J% = t{-^)~ never vanishes 



□2 Ak) 

"d^ ~ ''^ dX 

for t ^ 0. Hence, by the stationary phase expansion, is a lagrangian 
distribution whose principal symbol can be written as a smooth function 
of A : 

\Qr 

with Q = — jTA:(A)(^)~^. Since its amplitude vanishes precisely with x, we 
can deduce that the semiclassical wave- front of is (for fixed j (z 7^*) 

WFniZk) = {(A, t) G r2. ^ ^ -jrfc(A), x(A) / 0}. 

We still need to sum up all Zk{t;j,h) for j € Z*. For this we consider the 
localisation of Z. Without loss of generality, one can restrict to positive 
times. Let to > 0, e > 0, and let p £ CQ°{B{tQ,e)). There is no solution 
to (11) in the support of p for \j\ outside the interval 



hie) 



f tp-e tp + e 
Vsupj jrfcl' infj jr^l 



Making explicit the non-stationary phase argument, we can write, for any 
£ G N, 

Z,{t;j,h) = (^)7 e-^'^-^(^^+*^o'='(^^^))L^(a(c;n))dc, 
where L is the linear differential operator defined by 

(Lu)(c) = — 



j dc 

and a{-]K) € C^{I) admits an asymptotic expansion in non- negative powers 

of in the C°° topology. Let 6(c) = (1 + j Qc )~ Then h is uniformly 
bounded on / for \j\ > {tp + e)/ inf j |rfc|, and for any £ € N*, there exists 

^ Ce/j. 



a positive constant Ce, independent of j and h, such that 
Therefore, there exist constants Ci > such that 
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and we get, again when \ j\ > (io + e)/ infj \Tk 

\p{t)Zk{t;j,h)\^Ci 

Thus, for i^2, 

\p{t)Zi,{t;j,h)\^Ceh'. 



This shows that only a finite (independent of h) number of terms contribute 
to p{t)Z{t;h) modulo 0{h°°). Thus the (non)-stationary phase approxima- 
tions are jointly valid. Therefore Z{t; h) microlocally vanishes at any point 
that does not belong to C{I); this writes 

WFn{Z{-;h)) cC{I). 

More precisely, 

WFn{pZ{-,h)) C {(A,irfe(A)); A G /, U| € 4(e), A: = 1, . . . , iV} . 

Moreover, at a non-resonant point (A,jTfc(A)), no other period j'T^' can 
contribute, and Z{-;h) is a lagrangian distribution with principal symbol 
equal to Lf:{X;j, h). Since the set of resonant points is discrete, and WFji{Z) 
is closed in T*/, we must have WFfi{Z{-; h)) = C{I), which finishes the proof 
of the first step. 

Step 2. We are now left with a simple geometric inverse problem : given 
the set of periods how can one recover the number N of connected 

components ? 

Our strategy is to recover the fundamental periods |ri| , . . . , |TAr|. First 
of all, by Weyl's asymptotics (6), one obtains the a priori bound |rfc(A)| ^ 
7rp(X). Let R := maxjvrp. Then by assumption, the set of resonant points 
inside /x]0, R] is finite; therefore, one can always find a smaller, non-empty 
interval I C I such that there is no resonant point at all in /x]0, R\. 

We extract the periods r/j from Ci := C{I) fl (/x]0,i?]) inductively, as 
follows. 

1. Consider a point (Ai,ri) € £i with "minimal height" ti : V(A,r) G 

Cl,Ti ^ T. 

2. By the non-resonance assumption, the connected component of (Ai, ri) 
in Ci is the graph of a smooth function of the interval /. We denote 
this function by A i— > ti(A). 
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3. Consider the set 

£2 := A\{(A,iri(A)); A G /,j G Z*}. 

Again by the non-resonance assumption, L\ remains a union of non- 
intersecting smooth graphs. 

4. If L\ is empty, then = 1. Otherwise, start again by replacing by 
£1, and so on. If L]^ is empty, then = A; — 1. 

□ 

Remark 4.3 If we disregard symmetry issues, our assumption on the res- 
onant set is quite weak. For instance, one can allow the crossing of two 
periods to be flat (all derivatives are equal at a point A), simply because 
we put ourselves in a region with no crossing at all. However, it is easy 
to prove Step 2 with even weaker assumptions. For instance, it may work 
even if there are some open intervals of values of A which admits resonant 
pairs. It would be interesting to know whether Step 1 could hold in this case 
as well. It would then involve sub-principal terms in the Bohr-Sommerfeld 
expansion. A 

4.2 Singular fibres 

As we already mentioned, the following result comes for free. 

Theorem 4.4 Let Assumption A{P, J ^ I) hold, and let Eq £ I be a non- 
degenerate critical value of p. Assume also that p~^{Eq) contains only one 
critical point. Then from the knowledge of T,{P, J' , I) + 0{h'^) one can de- 
termine the topology of the singular foliation induced by p, in a saturated 
neighbourhood of p~^{Eq). 

Proof . Under these assumptions, the topology of the singular foliation 
induced by p in a saturated neighbourhood of p~^{Eq) is known to be com- 
pletely characterised by the type of the singularity [8, 24], which is deter- 
mined by Theorem 3.1. For the convenience of the reader, we briefly recall 
the two possible cases. 

1. The elliptic case. — The singular fibre p~^{Eq) is just a point and 
the foliation is homeomorphic to the one given by the Hamiltonian 
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2. The hyperbolic case. — The singular fibre is a circle with a transversal 
self- intersection (the figure eight). It separates a saturated neighbour- 
hood into three connected parts : two on one side, and one on the 
other side. It is homeomorphic to the foliation given by the Hamilto- 
nian H{x, ^) = .^^ + — x^, in a neighbourhood of i/~^(0). 

□ 

4.3 Global topology 

We say that a hamiltonian system p on the symplectic 2-manifold M is 
topologically equivalent to the hamiltonian system p on M if there is a 
homeomorphism if : M ^ M such that 

p = po ip. 

Notice that this implies that if respects the foliation, fibre by fibre. In 
particular, p and p have the same set of regular values and the same set 
of critical values. If / is an open interval, then two hamiltonian systems p 
and p are called topologically equivalent over / when they are topologically 
equivalent when restricted to the symplectic manifolds p~^{I), {p)~^{I). 

We call the topological type of a hamiltonian system the equivalence class 
of topologically equivalent systems. 

Theorem 4.5 Let Assumption A{P, J ^ I) hold, and assume that p has only 
nondegenerate critical values in some neighbourhood of I, such that any two 
critical points with the same singularity type cannot have the same image 
by p. Let ci < • • • < c„ 6e the critical values of p in L. Suppose that in 
each interval (cj, Cj+i), i = l,...,?i— 1, there exists a non-empty subinterval 
Li such that the set of resonant points in Li x M is discrete in M?. Then 
the knowledge of Ti{P, J , L) + 0{h?) determines the topological type of the 
hamiltonian system p over L . 

Proof. Upon a possible enlargement of /, one may assume that / = {Eq, Ei) 
for regular values Eq,Ei. Using symplectic cutting [13] or surgery [24], one 
may replace the phase space M'* by a compact symplectic manifold where 
p~^{L) is embedded. Then we apply the result of [8] that says that the 
topological type of p on M is determined by its Reeb graph : the set of leaves 
of the foliation, as a topological 1-complex. This graph is characterised by 
the relative positions of critical values, and the number of fibres between two 
consecutive critical values. The former is determined by the spectrum in / 
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thanks to Theorem 3.1, while the latter is determined for each i = 1, . . . , n— 1 
by the spectrum in /j, thanks to Theorem 4.2. This give the topological 
type of p, up to some homeomorphism of the Reeb graph itself. But since 
we know the precise values of p at singularities, we can in fact assume that 
this homeomorphism is the identity. □ 



5 Symplectic geometry 

The tools we've used so far give us the periods of the classical hamiltonian 
system, which is of course much more than a mere topological information. 
We show here that it is indeed sufficient to recover the full dynamics of the 
systems. 

We say that a hamiltonian system p on the symplectic 2-manifold M 
is symplectically equivalent to the hamiltonian system p on M if there is a 
smooth symplectomorphism ip : M ^ M such that 

p = po if. 

Thus, the dynamics of p on the levelset {p = E} is transported via ip to 
the dynamics of p on the levelset {p = E}. 

We call the symplectomorphism type of a hamiltonian system the equiva- 
lence class of symplectically equivalent systems. As before, one may restrict 
this equivalence to an interval / of values of p and p. 

Definition 5.1 Let (A,t) € / x M be a resonant point for p. Thus 

jTk{X) =j'Tk'{X) = -t 

for some j, j' ,k ^ k' . We say that this resonance is weakly transversal if 
there exists an integer n S N* such that the n-th derivatives of the periods 
are not equal : 

Theorem 5.2 Let Assumption A{P, J , L) hold, and suppose thatp has only 
nondegenerate critical values in some neighbourhood of I, such that any two 
critical points with the same singularity type cannot have the same image 
by p. Let ci < ■ ■ ■ < Cn be the critical values of p in L . Suppose that for 
each interval Ji := (q, Q+i), i = 1, . . . , n — 1, the set of resonant points in 
Ji X R is discrete. Finally assume that all such resonant points are weakly 
transversal. 

Then the knowledge ofT,[P,J',I) + 0{h?) determines the symplectic type 
of the hamiltonian system p over L . 
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Proof. We use the symplectic classification of [8, 17] using weighted Reeb 
graphs. Under our assumptions, the Reeb graph has vertices of degree 1 and 
3. A vertex of degree 1, a bout, corresponds to an elliptic critical value, while 
a vertex of degree 3, called a bifurcation point, corresponds to a hyperbolic 
critical value. At a bifurcation point we can distinguish one particular edge, 
called the trunk, corresponding to the side of the figure 8 with only one 
connected component. The two other edges are called the branches. A 
weighted Reeb graph is a Reeb graph each of whose edges is associated with 
a positive real number, its length, and such that each of the two branches of 
each bifurcation point is associated with a formal Taylor series {ie a sequence 
of real numbers). The hypothesis of the theorem allow for determining the 
topological Reeb graph via Theorem 4.5. Thus, the next step of the proof 
is to show how the numbers that constitute the weighted Reeb graph can 
be recovered from the spectrum. The final step is to obtain the symplectic 
equivalence in the sense that we have just defined above. 



The lengths. — Let Ck{Ji), for k = 1,...,N, be the connected com- 
ponents of p~^{Ji). Let Kk^i € C°°{Ck{Ji)) be an action variable for the 
regular lagrangian fibration P\c^{Ji)] it is unique up to a sign and an addi- 
tive constant. By definition the length of the edge corresponding to the set 
of leaves in Ck{Ji) is 



4, 



lim Kk,i{c) - lim Kk^i{c] 



(12) 



In a learned terminology, this is the Duistermaat-Heckman measure of Jj 
for the S^-action defined by K^ i, or, equivalently, it is the affine length of 
Ji endowed with its natural integral affine structure given by P\Cf,{Ji)- 

It follows from the local models for elliptic and hyperbolic singularities 
that this length is always finite. This is obvious at elliptic singularities, 
where the action has the form + At a hyperbolic singularity m, 
one can introduce a foliation function q such that, in some local symplectic 
coordinates around m, q = x^, and g > on the branches while g < on 
the trunk. Then the Duistermaat-Heckman measure has the form 

idfijiq) = {lnq + gj{q))dq on each branch (j = 1,2) ^^^^ 
1 (i/.f(g) = (2 In \q\ + g{q)) dq on the trunk, 

with some smooth functions g,gi,g2 satisfying 

Vp, 5^^)(0)=5?)(0)+5?^(0). 
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Under this form, the Taylor series of the functions g, gi, 52 at the origin are 
uniquely defined [17, 19]. 

Using the proof of Theorem 4.2, from the spectrum in I we can recover 
the periods rfc(A), A: = 1, . . . , A^, for A in any interval in Jj where the graphs 
of the periods don't cross. At a crossing the difficulty is to put the labels 
k correctly, so that the connected components Cfc(A) remain in the same 
Ck{Ji) when A varies. This can be overcome precisely thanks to the weak 
resonant assumption at each crossing, because each is C°° in Jj. This 
was the main issue. Now, fixing a point Aj € Jj, the action variable K^^i can 
be computed by the formula 



This gives the length of Ck{Ji) via equation (12). 

The Taylor series at the bifurcation points. — By definition, the 
sequences of numbers associated with a bifurcation point in the Reeb graph 
are the Taylor series of the functions gi, g2 (defined in equation (13)) at the 
origin. 

Let us show how to recover the Taylor series of g from the spectrum. 
The procedure is completely analogous for gi and g2- 

Thus, we consider a hyperbolic critical value Cj+i. We want to express 
the Duistermaat-Heckman measure on the trunk in terms of the principal 
symbol p. By a theorem of Colin de Verdiere and Vey [4], there exist local 
symplectic coordinates (x, at the hyperbolic point, and a smooth, locally 
invertible function / : (M, Cj+i) (1^,0) such that 



For notational purposes, one may assume that /'(cj+i) > 0, which amounts 
to say that the trunk is sent by p to A < Cj+i. Then from (13), for A close 
to Ci+i, A < Ci+i, 



On the other hand if the connected component corresponding to the trunk is 
Ck{Ji), one has by definition of the Duistermaat-Heckman measure dii{X) = 
Tk{X)d\. Therefore 




f{p) =x^ = q. 



d/i(A) = (21n|/(A)|+go/(A))/'(A)dA. 



Tfc(A) = /'(A) (2 In |/(A) \+go /(A)) = 2/'(A) In | A - q+i | + h{X) 
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for some smooth function h at X = Q+i. There, using Taylor's formula, we 
have written /(A) = a(A — Cj+i) + (A — Cj+i)^/(A), with a > and / smooth 
at Cj+i, and hence 



hiX) = 2 fix) In a + (A - Q+i)/(A) + fiX)g o /(A). 



(14) 



This shows that h is smooth for A close to Cj+i. 

It is easy to see that any smooth function i;^ in a neighbourhood of the 
origin such that (p{t) In t extends to a smooth function at t = must be flat. 
Hence the knowledge of Tk{X) for A < Cj+i completely determines the Taylor 
series of /'(A) (and hence /(A)) at A = Cj+i. 

Then one can recover the Taylor series of h using 



Finally, from (14) and the fact that / is locally invertible, one can recover 
the Taylor series of g at the origin. 

Symplectic equivalence. — We have proven that the weighted Reeb 
graph is determined by the spectrum. By Toulet's classification [8, 17], if 
two such systems {M,p) and {M,p) have the same weighted Reeb graph, 
there exists a symplectomorphism ip : M ^ M such that p and po (p define 
the same singular foliation on M (ip indices a homeomorphism of the leaf 
space, fixing the vertices). If we assume that the operators P and P have 
the same spectrum (modulo h"^) and fulfil the requirements of the theorem, 
then we also know that p and p o (p share the same set of critical values 
Cj. The fact that p and p o ip define the same foliation implies that for 
each connected component Ck{Ji), there exists a smooth, invertible function 
f : Ji ^ Ji such that 



Since the singular fibres at the ends of Ck{Ji) are fixed by ip, f must be 
increasing, and thus extends to a homeomorphism of Jj. 

As we already saw, the spectrum also determines the periods at a given 
energy E = X. Hence for X G Ji, rfc(A) = ffc(A). Since is integrable at 
Cj+i, we can define action integrals for A < Cj+i as : 



MA) =rfc(A)-2/'(A)ln|A-Q+i 



VA < Q+i 



p = f opoip 



on Ck{Ji). 



(15) 
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We have Kk^i{\) = Kk^i{\). On the other hand, the action is a symplectic 
invariant of the fohation. Prom (15) on can compute the action on the curve 
^{Ck{f{m = Ck{\) : i^fc,i(/(A)) = i^fc,i(A) + const. Therefore 

Since Tk does not vanish in Jj, K/^^i is strictly monotonous on Jj. Therefore 

/(A) = A, VA e Ji. 
Thus p = po (f on each Ck , and by continuity 

p = po if on M. 

This finishes the proof of the theorem. □ 
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